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We describe a simple model for the melting and optical properties of a DN A/gold nanoparticle 
aggregate. The optical properties at fixed wavelength change dramatically at the melting transition, 
which is found to be higher and narrower in temperature for larger particles, and much sharper than 
that of an isolated DNA link. All these features are in agreement with available experiments. The 
aggregate is modeled as a cluster of gold nanoparticles on a periodic lattice connected by DNA 
bonds, and the extinction coefficient is computed using the discrete dipole approximation. Melting 
takes place as an increasing number of these bonds break with increasing temperature. The melting 
temperature corresponds approximately to the bond percolation threshold. 



PACS numbers: 61.43.Hv, 78.67.-n, 82.60.Qr, 87.15.-v 



The optical properties of metallic nanoparticles have 
been investigated intensively over the last decades [1,2]. 
Recently, this investigation has expanded to include 
so-called functional metallic nanoparticles [3]. Among 
these, there is a particular interest in DNA-modified gold 
nanoparticles. Such particles can form complex aggre- 
gates which may melt at high temperatures. The parti- 
cles and their aggregates may have a variety of applica- 
tions, e. g., in biological detection, which may be possible 
by using the optical and electrical sensitivity of the aggre- 
gates [4-6]. Numerical model calculations of the optical 
properties of DNA modified gold nanoparticle aggregates 
show general agreement with experiments, including such 
features as the following: (i) for isolated gold nanopar- 
ticles in suspension, there is a strong "surface plasmon" 
absorption in the visible; and (ii) this absorption max- 
imum broadens and red-shifts when the cluster radius 
becomes comparable to the wavelength [7,8]. However, 
some important physical details have not been explained. 
For example, in the gold/DNA nanoparticle system, the 
melting transition of a gold nanoparticle aggregate has a 
much narrower temperature width and occurs at a high 
temperature than that of a single DNA link [5,9]. Also, 
the dependence of aggregate melting temperature on par- 
ticle size has not yet been explained [10]. 

In this Letter, we model this novel "melting" transi- 
tion of a gold/DNA nanoparticle aggregate. Our model 
accounts for most experimentally observed features, in- 
cluding (i) the small temperature width of the aggregate 
melting transition, compared to that of a single DNA 
link; (ii) the particle-size-dependence of the melting tem- 
perature T m ; and (iii) the temperature-dependence of the 
optical extinction coefficient C ext (X,T) at wavelength A 
and temperature T. 

We first describe our model for the low-temperature 
morphology of the aggregate. In the simplest version of 
our model, the low-T cluster is taken simply as a col- 
lection of identical gold nanoparticles (each of radius a) 



which fill the sites of a simple cubic lattice of lattice con- 
stant d (d > 2a). The cluster is assumed to be a cube 
of edge L, containing N pax = {L/d) 3 gold nanoparticles. 
We have also investigated the melting and the optical 
properties assuming that the low-T cluster is a fractal 
aggregate. 

To describe the melting of this aggregate, we assume 
that each nanoparticle carries exactly N s single DNA 
strands. (This choice still retains the essential features of 
the melting.) The cluster exists at low T because a chem- 
ical reaction converts two single strands S into a double 
strand D. In this Letter, we adopt a simplified two-state 
model for this reaction [11,12], described by the relation 



S + S^D. 



(1) 



For short DNA strands (12-14 base pairs), this model 
describes melting well [13]. Each gold nanoparticle has 
z nearest neighbors (z — 6 for a simple cubic lat- 
tice). Hence, there are N paT z/2 bonds joining adja- 
cent nanoparticles. We assume that a (temperature- 
dependent) fraction p(T) of the single DNA strands form 
double strands by the reaction (1). In order to calculate 
the fraction p c f[(T) of bonds which contain at least one 
double strand, we adopt the following model. First, we 
assume that exactly N s /z of the single strands on a given 
nanoparticle are available to bond with any one of its z 
nearest neighbors. Thus, we assume that the maximum 
number of links that could be formed between any two 
particles is N s /z. The probability that no link is formed 
is then taken to be 



1 - PcB (T) = [1 - p(T)} 



N„/z 



(2) 



The criterion for the melting temperature T c is that 
Pes(T c ) = p c , p c being the bond percolation threshold 
[14] for the lattice considered, at which an infinite con- 
nected path of double DNA strands first forms. (For 
example, p c ~ 0.25 on a very large simple cubic lattice.) 
Thus, at percolation, 



1 



[l-p(T c )] N -/' = l- 



Pc 



(3) 



This equation implicitly determines T c in terms of z, p c , 
and N s [15]. 

p(T) itself is determined by a chemical equilib- 
rium between the single-strand and double-strand DNA 
molecules, which are attached to the gold nanoparticlcs. 
The chemical equilibrium condition corresponding to (1) 
is 



P(T) 



(4) 



where K{T) is a suitable chemical equilibrium constant, 
and Ct is the molar concentration of single DNA strands 
in the sample (in our model all are attached to the 
nanoparticles). Since < p(T) < 1, the physical solu- 
tion to cq. (4) is p(T) = 1 + \ (K' - Vif'- 2 + AK>) . 

Note that, according to eq. (3), p(T c ) decreases with in- 
creasing 7V S and, hence, with increasing particle radius. 
Since p(T) decreases monotonically with T, T c should 
thus be an increasing function of a, as reported in exper- 
iments [10]. 

In Fig. 1, we plot p e ff(T) for several a. We assume 
N s oc a 2 , set z = 6 and use the experimental re- 
sult that N s = 160 when a = 8 nm [16]. We have 
also assumed the simple van't Hoff behavior K(T) = 
exp[— AG/fcsT], with a Gibbs free energy of formation 
AG(T) = ci(T - T M ) + c 3 (T- T M ) 3 , choosing the values 
of Ci, C3, and Tm to be consistent with experiments on 
these DNA molecules. 

To calculate the T-dependent optical properties, we 
have considered the melting of two slightly different low- 
temperature aggregates. In the first, we assume that, at 
low T, the aggregate consists of a simple cubic collection 
of iV par = (L/d) 3 gold nanoparticlcs on a simple cubic 
lattice, as described above, with all bonds occupied by 
DNA double-strand links. To generate a specific sam- 
ple with a given p c fj (T), we randomly remove links with 
probability 1 — p c tf(T), then identify the separate clus- 
ters, using a simple computer algorithm [14]. We neglect 
gravitational forces, which may be important in some 
experimental circumstances. If the aggregate consists of 
two or more clusters, we simply place these aggregates in 
random positions and orientations within a larger bound- 
ing box (usually of edge lOOd), taking care that the in- 
dividual clusters do not overlap. The resulting geometry 
is shown schematically in Fig. 2. We have also carried 
out the same procedure to simulate the melting of a sam- 
ple formed by reaction-limited cluster-cluster aggregation 
(RLCA) [17]. A typical RLCA cluster is shown in Fig. 
2 (d), and represents a possible fractal aggregate which 
might be produced by certain random growth processes 
at low T [18]. 

We calculate the optical properties of this sample us- 
ing the Discrete Dipole Approximation (DDA) [19,20]. 



The sample is modeled as a collection of separate ag- 
gregates whose extinction coefficients are computed indi- 
vidually, then added. Each aggregate consists of many 
identical nanoparticles, which have complex frequency- 
dependent dielectric constant e(u>), and polarizability 
a(u) related to e(w) by the Clausius-Mossotti equation 
[21] e(w) = 1 + 47rna/[l - (47rn/3)a], where n = {d/Lf. 
The resulting expressions for the induced dipole moment 
Pi of the i th sphere, and the corresponding expression for 
the extinction coefficient C ex t(k) at wave number k, are 
given in Ref. [7]. 

In our case, each cluster consists of a number of DNA- 
linked individual gold nanoparticles. In our calculations, 
we do not include the optical properties of the DNA 
molecules, since these absorb primarily in the ultravi- 
olet [8]. We use tabulated values of the gold complex 
index of refraction [22,23], then calculate C ex t for each 
cluster using the DDA. To improve the statistics, we av- 
erage C ex t for each cluster over possible orientations. We 
then sum the averaged extinction coefficients of all the in- 
dividual clusters to get the total extinction coefficient of 
the suspension. This method is justified when the sus- 
pension is dilute. 

In Fig. 3, we show C ext (X,T) versus A for the gold- 
DNA cluster at several values of p c ff(T), assuming a = 20 
nm. p c ff = represents a dispersion of iV par = {L/d) 3 in- 
dividual nanoparticles, while p c ff = 1 represents a simple 
cubic lattice of connected nanoparticlcs. The calculated 
C ex t(X,T) at each A changes strikingly for p e gf(T) <~ p c . 
For small p e ff, there is a clear extinction peak near 520 
nm. This peak corresponds to the wavelength of the sur- 
face plasmon resonance (SPR) in individual Au nanopar- 
ticles. As Peff increase, this peak first red-shifts, then 
greatly broadens, as the aggregate melts. The calculated 
peak is, however, shifted much more and is much broader 
than experiment. However, as shown in inset, if we as- 
sume that the low-T aggregate is an RLCA fractal, the 
calculated peak shift is consistent with experiment [18]. 

In Fig. 4, we show C ex t(X,T) at fixed X = 520 nm 
(close to the isolated-particle SPR), versus T, for sev- 
eral particle sizes, assuming that the low-T aggregate is 
a simple cubic cluster with N — 1000 particles, as shown 
in Fig. 2. For each a, the extinction increases sharply at 
a characteristic T, corresponding to the melting of the 
aggregate for that a; at this T, the absorption due to the 
SPR increases sharply. 

In the inset of Fig. 4, we compare C ex t(X, T) for a reg- 
ular and an RLCA cluster of particles of 20 nm radius 
at A = 520 nm, both for N = 1000 particles. Although 
the RLCA cluster has a slightly broader melting transi- 
tion, as manifested in C ex t(X,T), than does the regular 
lattice, both sets of data show a much sharper melting 
transition than that of a single DNA link. Also, although 
our normalized C ext (X, T) is calculated for the aggregates 
at 520 nm, we expect similar behavior at 260 nm. (We 
have not carried out calculations at this A mainly because 
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we have not included the DNA absorption properties.) In 
any case, the experimental melting curves at 260 and 520 
nm are very similar [10]. 

Our calculated extinction coefficients are strikingly 
similar to recent experimental results [5,8,10]. In par- 
ticular, both experiment and calculation give a sharp 
increase in C ext (X, T) at fixed A, as T increases past a 
critical temperature, which we interpret as the melting 
temperature T m . We also find, in agreement with exper- 
iment [5,8], that melting occurs over a much narrower 
range of T in the aggregate than for a single bond, and 
that the melting occurs at higher T for larger particles 
[10]. The crucial point is that p e &(T) is a much sharper 
function of T than is p(T), and this feature would not 
be affected by slight changes in the model (such as con- 
sidering a body-centered-cubic rather than a simple cu- 
bic cluster). In the present model, the melting transi- 
tion is in the universality class of bond percolation. This 
universality class might change because of the constant 
dissolving and reforming of the bonds of the aggregate 
(a possibility omitted from our model). This possibility 
should be studied further theoretically and experimen- 
tally. However, the optical properties may well be robust, 
in the sense that C ex t(T) would be little affected even if 
this feature were included, since the key ingredient is the 
multiple links per bond as described above. 

In summary, we have developed a simple model for the 
melting of aggregates of gold nanoparticles and DNA, 
and have calculated the T-dependent optical properties of 
these melting aggregates, using the DDA. We find that, 
at fixed A, melting is accompanied by dramatic changes 
in the extinction coefficient C ext (X). These calculated 
changes occur over a much narrower temperature range 
than that over which interparticle links themselves melt, 
the temperature width is smaller for particles with larger 
radius, and the melting occurs at higher temperatures 
for the larger nanoparticles than for the smaller ones. 
All these effects are in good agreement with experiment 
[5,8,10]. It would be of great value if the predicted In- 
dependent structure of the clusters could be probed ex- 
perimentally, e. g. by static light scattering, as has been 
done in other contexts [24] . Finally, we mention that the 
present model can be applied to other similar nanocom- 
posites [25]. 
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FIG. 2. Schematic of the melting of a gold-DNA cluster, 
for two different models discussed in the present paper. In 
the first model, (a) at low T the cluster is described as an 
L x L x L simple cubic aggregate of lattice constant d [here 
L/d — 10]. (b) As the temperature T increases, some of the 
bonds break. The fraction of bonds present is p c s(T). In 
(b), p cS = 0.5 > p c (L), where p c is the weakly //-dependent 
percolation threshold. A percolation cluster of linear dimen- 
sion ~ L still exists. At a still higher T [shown in (c)], 
Pes(T) = 0.2 < p c (L), and only small clusters remain. In 
our calculation, the clusters are positioned and oriented at 
random in the bounding box, taken as a cube of edge lOOd. 
(d) Alternate model for low-T sample (p B s(T) — 1): fractal 
cluster formed by reaction-limited cluster-cluster aggregation 
[RLCA], with fractal dimension df ~ 2.1. 
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FIG. 1. Plot of p e g(T) versus T for several different choices 
of particle radius a, as indicated. Also plotted is p(T), the 
probability that a given DNA strand is part of a double strand 
at T. 
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FIG. 3. Calculated extinction coefficient C ex t(^,T) as a 
function of wavelength A, plotted for several values of p e g, 
particle radius 20 nm, as indicated in the legend. p e s = 1 cor- 
responds to 1000 gold nanoparticles on a simple cubic lattice 
of edge L = 480 nm. The sample is that shown schematically 
in Fig. 2 (a) - (c). Inset: C ex t(\,T) for the RLCA aggregate 
at p e ff = 1 (curve A), as well as for simple cubic aggregate 
at p c g — 1 (curve B), and for individual gold nanoparticles 
Pes = 0. (curve C). 
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FIG. 4. Normalized extinction coefficient C ext {\,T) for 
A = 520 nm, plotted versus T for several particle radii o, 
assuming that the low-T aggregate is an N = 1000 simple 
cubic cluster, as shown in Fig. 2. Inset: Normalized extinc- 
tion coefficient C ex t(A,T), versus T for A = 520 nm, plotted 
for a 1000-particle gold/DNA aggregate assuming that the 
low-temperature sample is a simple cubic cluster (open cir- 
cles) or an RLCA cluster (open squares). The solid curve is 
a plot of 1 — p(T) for a single DNA duplex with the same 
concentration Ct as the above two curves. The dotted curve 
represents 1 —p{T) for a single DNA duplex but with a much 
higher Ct than for the other curves of the inset. The plots 
of 1 — p(T) for a single link closely resemble C ex t(A,T) as 
measured at A = 260 nm for a single duplex[ll]. 
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